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CONVOLUTIONAL  CODING  OPTIONS  FOR  MFSK/FH 
SIGNALING  ON  RALEIGH  FADING  AND 
PARTIAL  BAND  INTERFERENCE  CHANNELS 

1.  INTRODUCTION 


It  is  well-known  that  forward  error  correction  coding 
offers  significant  performance  advantages  when  applied  to  the 
reception  of  signals  in  coherent  additive  white  Gaussian  noise 
( AWGN)  channels.  Por  example,  with  a  constraint  length  7 ,  rate 
1/2,  convolutional-code  with  soft  decision  decoding,  the  E./N 
required  for  a  10~  bit  error  probability  is  reduced  by  about°5 
dB  relative  to  the  uncoded  system.  In  1975,  Viterbi  and  Jacobs 
[1]  showed  that  much  larger  coding  gains  are  achievable  on 
noncoherent  channels  with  either  amplitude  fading  or  additive 
partial  band  Gaussian  interference.  In  fact,  even  with  modest 
coding  protection,  coding  gains  greater  than  30  dB  are  possible 
for  fading  or  worst  case  partial  band  Interference  channels. 

In  this  report  we  investigate  several  convolutional  coding 
options  which  are  applicable  to  the  noncoherent  fading  or 
partial  band  Gaussian  channel.  These  results  are  useful,  for 
example,  on  an  HP  channel  where  multiple  ionospheric  reflections 
lead  to  a  loss  of  phase  information  as  well  as  introducing  a 
fading  amplitude  characteristic.  In  other  instances,  channels 
exist  with  partial  band  interference  caused  either  by 
unregulated  other-users  or  by  intentional  jammers. 

The  modulation  technique  considered  here  is  multiple 
frequency  shift  keying  with  frequency  hopping  (MPSK/PH) .  With 
orthogonal  tone  spacings,  MFSK  is  an  optimal  M-ary  modulation 
for  noncoherent  channels.  Frequency  hopping  and  (possibly) 
interleaving  are  assumed  to  be  used  to  make  the  channel 
memoryless,  that  is,  to  produce  statistical  independence  from 
symbol  to  symbol  in  a  sequence  of  received  symbols.  In  the 
following  section  the  uncoded  performance  of  this  modulation 
technique  is  given  for  the  two  channel  models  in  question. 

2.  PERFORMANCE  OF  UNCODED  MFSK/FH  ON  RAYLEIGH  FADING  AND 
PARTIAL  BAND  GAUSSIAN  CHANNELS 

For  orthogonal  MFSK  the  M-ary  signal  alphabet  is  chosen  as 
the  set  of  tones: 


s1  (t) 


COs(2irfjt  +  9 ),  0stsT$ 


(1) 


where 

Es  is  the  received  symbol  energy, 

Tg  is  the  symbol  duration, 

and  es  is  a  uniformly  distributed  random  phase. 
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The  frequencies  are  chosen  to  insure  mutual  orthogonality  of  the 
signals,  so  that  the  tone  spacings  are 


^,1*1,2,  ...  M-l 


(2) 


For  a  symbol  to  carry  k  bits  of  information, 
follows  that 


E 


b 


we  use  M=2 


k 


and  it 


(3) 


where  is  the  received  energy  per  bit. 

For  a  Rayleigh  fading  channel  when  signal  s.(t)  is  sent,  the 
received  signal  is  1 


r(t)  ■  r 


cos(2irf..t+0)  +  n(t),  0£t*Ts 


(4) 


where  n(t)  is  an  AWGN  process  with  two-sided  power  spectral 
density  Watts  per  Hertz,  and  r  is  a  Rayleigh  random  variable 
with  normalized  probability  density  function 


p(r)  »  2re“r  ,  r>0  (5) 


2 

With  this  normalization,  r  E  is  the  received  energy  random 
variable  whose  average  value8 is  Eg. 


The  uncoded  performance  of  MFSK  in  a  memoryless  Rayleigh 
fading  channel  is  well-known  [2].  The  symbol  error  probability 
Pg  is  given  by 


P 


s 


1 


(n+1 )+n  m- 
o 


(6) 


Furthermore,  with  orthogonal  signaling,  bit  and  symbol  error 
probabilities  are  related  by 


pb  ’  ZftPTT  % 


(7) 


Thus  from  (6) ,  (7) ,  and  (3)  we  have 

V ■  rfTT  g  (v)  (")  —‘-5 

^  '  '  '  '  (n+1  )+nk  ^ 

The  curves  representing  (8)  are  plotted  in  Figure  1  for 
M®2,4,...32.  We  see  that  for  the  case  of  M*2,  an  tb/No  of 
approximately  50  d8  is  required  to  achieve  a  bit  error 
probability  of  10“5 .  Increasing  M  offers  little  improvement;  in 
fact  it  has  been  shown  [2]  that  as  M—  oo ,  the  improvement  over 
the  binary  case  is  only  2/ln2  (4.6  dB) . 

For  uncoded  MFSK/FH  on  a  partial  band  Gaussian  interference 
channels  the  result  is  similarly  poor.  The  partial  band  channel 
is  characterized  by  constant  density  AWGN  over  a  fraction  of  the 
total  hoppinq  transmission  bandwidth.  Thus,  the  noise  spectrum 
has  density  N 0/p  over  a  fraction  p  of  the  band  (where  0<p<l)  and 
is  zero  elsewhere  (oyer  a  fraction  1-p  of  the  band) .  We  assume 
that  the  M  candidate  tone  slots  for  each  symbol  are  either  all 
interfered  with  or  they  are  all  noise  free.  As  a  worst-case 
condition,  we  consider  only  the  situation  where  the  parameter  p 
is  chosen  so  as  to  maximize  the  resulting  probability  of  error. 

For  uncoded  MFSK/FH  in  a  worst  case  partial  band  Gaussian 
channel  the  results  (over  the  range  of  interest)  have  been  found 
in  [3]  as 


P 


b 


b 

VN0 


(9) 


where  the  numerator  b  is  a  constant  depending  on  the  parameter 
M,  given  as  follows: 

M  b 

1  .17 

4  .23 

8  .20 

16  .18 

32  .17 

The  curves  representing  (9)  are  plotted  in  Figure  2.  The 
results  are  about  4  dB  better  than  the  corresponding  Rayleigh 
fading  curves,  but  the  performance  is  again  poor  because  the 
dependence  of  Pb  on  Eb/No  is  inverse  linear  (rather  than 
exponential  as  in  the  AWGN  channel) .  For  each  of  these  two 
channels  we  shall  see  that  coding  can  provide  a  substantial 
performance  improvement. 
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.  CODES  AND  THEIR  WEIGHT  STRUCTURE 

In  this  report  we  consider  the  performance  of  various 
convolutional  codes  on  M-ary  orthogonal  channels  with  Rayleigh 
fading  or  partial  band  Gaussian  interference.  The  general  block 
diagram  for  this  system  is  shown  in  Figure  3  .  We  restrict  our 
consideration  to  convolutional  codes,  and  assume  that  soft 
decision  Viterbi  decoding  is  used  throughout. 


CODING  CHANNEL 


Fig.  3  —  General  modulation/coding  configuration 


The  specific  case  of  rate-1/2  codes  with  4-ary  orthogonal 
modulation  shall  be  considered  in  detail.  The  techniques  used, 
however,  are  applicable  to  all  M-ary  modulations  and  for  various 
code  rates.  The  codes  considered  fall  into  two  groups:  binary 
convolutional  codes  (that  accept  single  bit  inputs),  and 
quaternary  convolutional  codes  (that  accept  pairs  of  bits  at  the 
input).  The  former  are  designated  by  their  constraint  length  K 
and  we  shall  consider  the  cases  K*3,4,5,6,  and  7.  From  the 
latter  class  we  shall  investigate  the  dual-2  and  triple-2  codes, 
which  have  equivalent  binary  constraint  lengths  K**4  and  K*6 
respectively. 

These  seven  encoders  are  shown  in  Figures  4  through  10.  All 
tap  connections  (code  generators)  are  optimized  so  that  the 
codes  will  have  maximum  (or  nearly  maximum)  distance  properties 
with  respect  to  4-ary  channels  [4,5]. 
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In  evaluating  performance  of  a  convolutional  code,  the  only 
code  property  of' importance  is  its  weight  structure.  The  code 
weight  structure  gives  the  number  of  decoded  errors  associated 
with  erroneous  decoding  paths  of  various  distances,  starting 
with  the  code's  minimum  distance.  Because  of  the  linear  (group) 
property  of  convolutional  encoders,  all  error  calculations  may 
be  based  on  the  assumption  that  the  code  has  an  input  which  is 
the  all-zero  binary  sequence. 

In  computing  distance  it  must  be  kept  in  mind  that  we  are 
dealing  with  an  M-ary  orthogonal  (equidistant)  signal  set.  A 
4-ary  modulator  performs  a  mapping  of  code  bit  pairs  into 
channel  signals  as  follows: 


00 

01 

10 

11- 


Since  S] ,  S2,  and  s 3  are  equidistant  from  so  in  the  Euclidean 
signal  space  of  the  channel,  we  see  that  these  signals  must  be 
treated  as  equal  adversaries  to  the  correct  signal  sg  at  the 
demodulator  output.  Distance  from  the  all-Sg  sequence  of 
signals  for  this  channel  is  the  number  of  non-sg  signals 
produced  by  the  modulator  as  it  maps  a  sequence  of  encoder 
output  code  bits  to  a  sequence  of  modulator  output  signals.  The 
total  number  of  l's  at  the  encoder  input  which  produce  all  of 
the  sequences  of  distance  d  is  called  the  weight  of  the  code  at 
distance  d.  A  list  of  the  weights,  starting  at  the  shortest 
(minimum)  distance  is  called  the  weight  structure  of  the  code. 

The  search  procedure  for  the  code  weight  structure  is  described 
in  Appendix  A.  A  table  of  the  code  weight  structure  for  the 
four  shortest  distances  for  each  code  considered  is  given  in 
Table  I. 
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Triple-2 


K*7 


K*6 


K=5 


d 

Hd 

d 

**d 

d 

«d  * 

Sd 

7 

7 

6 

9 

6 

12  5 

3 

8 

39 

7 

14 

7 

40  6 

15 

9 

104 

8 

62 

8 

144  7 

22 

10 

352 

9 

212 

9 

488  8 

196 

K-4 

Dual 

2jL_ 

K-3 

d 

Hd 

d 

**d 

4 

1 

4 

4 

3  1 

5 

14 

5 

16 

4  4 

6 

21 

6 

56 

5  12 

7 

44 

7 

176 

6  56 

Table 

I. 

Weight 

4 

Structure 

(First 

Four  Entries) 

For  Codes  On  4-ary  Channels. 
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4.  PERFORMANCE  BOUNDS 


Upper  bounds  on  the  error  probability  performance  of 
convolutional  codes  with  MFSK  systems  can  be  found  by  using  the 
union-Chernof f  bound  technique  [6],  Results  obtained  from  this 
method  are  usually  about  1  dB  higher  than  the  exact  result  for 
bit  error  probabilities  less  than  10**3.  in  the  union-Chernof f 
bound  technique,  we  consider  all  candidate  error  paths  remerging 
at  a  node  on  the  correct  (all-zero)  path  of  the  decoding  trellis 
(state  vs.  time  diagram) .  This  is  typically  illustrated  in 
Figure  11. 


time 


state 


Fig.  1 1  —  Paths  on  a  decoding  trellis 


A  union  bound  is  employed  because  the  events  corresponding 
to  choosing  the  various  incorrect  paths  are  not  mutually 
exclusive.  The  calculation  is  facilitated  if  we  overbound  the 
probability  of  deciding  on  an  incorrect  path  by 

P(A1+A2+A3+.  .  .)  <  PfAjJ+PfAjJ+PfAj)*.  .  .  (10) 

where  A, ,  Alt  A3  ...  are  the  error  events  corresponding  to 
individual  incorrect  paths. 

The  importance  of  the  Chernoff  bound  lies  in  the  fact  that 
the  Chernoff  bound  is  itself  a  moment  generating  function.  (See 
Appendices  B  and  C.)  If  the  probability  of  error  for  a  single 
use  of  a  memoryless  channel  has  a  Chernoff  bound  B,  then  the 
Chernoff  bound  for  n  uses  of  the  channel  is  Bn .  In  the  decoding 
trellis,  the  probability  of  choosing  an  incorrect  path  of 
distance  d  is  therefore  Chernoff  bounded  by  B 4.  All  error  paths 
of  distance  d  have  a  collective  total  of  N.  ones  in  the  decoded 
bit  sequence,  where  N.  is  the  code  weight  as  defined  in  the 
preceding  section.  Tne  union-Chernof f  bound  on  the  decoded  bit 
error  probability  is  given  by 

H,  ~  ir  T  Mi  BA 

drrtin 

for  codes  with  binary  inputs  (K  ■  3  through  K  *  7)  and  by 


CO 


(11-b) 

for  codes  with  quaternary  inputs  (dual-2  and  triple-2) .  In 
equations  (11),  the  factor  of  1/2  in  front  of  the  summation  is 
an  additional  tightening  factor  which  can  be  applied  to  Chernoff 
bounds  under  fairly  general  conditions  [7).  In  equation  (11-b), 
the  code  weight  Nd  is  divided  by  a  normalizing  factor  of  2  to 
take  into  accountathe  fact  that  inputs  to  the  dual-2  and 
triple-2  code  are  pairs  of  bits  [8]. 

Equations  (11)  have  been  plotted  in  Figs.  (12)  and  (13) 
using  the  Chernoff  bound  B  derived  in  Appendices  B  and  C  for  the 
Rayleigh  fading  and  partial  band  Gaussian  channels.  The  values 
of  N.  are  those  given  in  the  preceding  section.  All 
calculations  are  truncated  after  the  first  four  terms. 
Furthermore,  in  the  partial  band  case,  it  is  assumed  that  all 
symbols  on  an  error  path  must  be  jammed  in  order  to  produce  an 
error.  This  is  the  "jammer  state  known”  assumption,  and  it  is 
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A 


frequently  used  in  analyses  [1]. 


By  comparing  Figures  (12)  and  (13)  with  Figures  (1)  and  (2) 
we  segethat  the  coding  improvement  (with  soft  decision  coding) 
at  10~3  error  probability  ranges  between  27  dB  and  35  dB  from 
the  weakest  to  the  strongest  codes.  In  the  next  section  we 
consider  the  decoder  complexity  for  the  codes  considered. 


5 .  DECODER  COMPLEXITY 


The  most  important  complexity  consideration  in  a  maximum 
likelihood  decoder  (employing  the  Viterbi  algorithm)  is  the 
number  of  pairwise  comparisons  required  in  the  decoding 
algorithm.  This  dictates  the  maximum  speed  for  serial 
processing,  or  computational  complexity  for  parallel  processing. 

For  a  binary  convolutional  encoder  of  constraint  length  K, 
the  number  of  encoder  states  is  2^-'  and  at  each  state  a  single 
comparison  is  required  for  each  input  bit.  For  dual-2  and 
triple-2  encoders  of  equivalent  (binary)  constraint  length  K»4 
and  K*6  respectively,  the  number  of  states  is  2  K-2,  but  at  each 
state  three  comparisons  are  required  to  determine  the  largest  of 
four  competing  path  metrics.  Since  the  encoder  accepts  two-bit 
inputs,  the  number  of  comparisons  at  each  node  per  input  bit  is 
3/2. 


In  Table  II  the  number  of  comparisons  per  bit  is  presented 
for  each  code  considered  in  this  report.  Using  this  as  a  simple 
complexity  measure,  we  see  that  the  complexity  of  the  decoder 
increases  as  the  performance  of  the  coding  system  improves.  (The 
single  exception  to  this  is  the  triple-2  code  which  slightly  outperforms  the 
K=6  code,  although  having  fewer  comparisons  per  bit). 


Code 

#  states 

fcomp. /bit/state 

u3*M‘iU*!u 

K-7 

64 

1 

64 

K*6 

32 

1 

32 

Triple-2 

16 

3/2 

24 

K=*5 

16 

1 

15 

K-4 

8 

1 

8 

Dual-2 

4 

3/2 

6 

K-3 

4 

1 

4 

Table  II. 

Decoder  Complexity  Measure  (comparisons 

per  bit) 

i 
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RESULTS  AND  CONCLUSIONS 


In  this  report  we  have  demonstrated  the  importance  of 
coding  in  transmission  of  MFSK/FH  signals  on  Rayleigh  fading  or 
partial  band  Gaussian  interference  channels.  For  a  probability 
of  error  of  10"'  the  coding  gain  of  a  K=7,  rate-1/2  code  on 
either  of  these  channels  is  approximately  35  dB.  For  the 
weakest  code  considered  (K*3) ,  the  coding  gain  is  nearly  27  dB. 

Using  a  simple  measure  of  decoder  complexity  (comparisons 
per  bit)  it  was  found  that  complexity  may  be  reduced  with  a 
resulting  loss  in  performance.  One  case  of  special  interest  is 
the  triple-2  code  which  has  performance  (at  10-5)  which  is 
within  1.5  dB  of  the  best  code  considered  (K*7  binary)  but  has 
complexity  which  is  less  than  one  half  of  that  of  the  K=7  code. 

Finally,  it  might  be  pointed  out  that  the  performance 
results  for  combined  Rayleigh  fading,  worst  case  partial  band 
Gaussian  noise  channels  are  implicit  in  the  results  given.  It 
has  been  shown  [9]  that  the  worst  case  Gaussian  jammer  for  a 
Rayleigh  fading  channel  is  a  broadband  jammer,  so  that  the 
Rayleigh  fading  results  are  also  applicable  to  the  combined 
case. 
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APPENDIX  A.  COMPUTATION  OF  CODE  WEIGHT  STRUCTURE 


The  weight  structure  of  a  convolutional  code  is  readily 
available  from  the  generating  function  for  any  given  code. 
However,  direct  calculation  of  the  generating  function  becomes 
difficult  for  codes  having  constraint  lengths  of  five  or  more, 
since  the  number  of  states  grows  exponentially  with  the 
constraint  length.  An  alternate  approach  is  to  implement  the 
encoder  and  calculate  the  weights  associated  with  the  various 
distance  error  paths.  The  path  errors  are  characterized  by  the 
codeword  distances  and  the  weights  are  determined  from  the 
number  of  input  "ones”  associated  with  each  codeword  distance. 
This  approach  is  equivalent  to  evaluating  the  generating 
function  at  some  selected  finite  sequence  length. 

A  computer  program  was  developed  to  compute  the  code 
weights.  It  does  so  by  implementing  the  coder  with  a  K-stage 
shift  register  and  2n  mod-2  adders  where  n  is  the  number  of 
4-ary  output  symbols  of  the  encoder.  The  shift  register  is 
shifted  b.  *  l(or  2)  bits  at  a  time  to  allow  binary  (or 
quaternary)  input  data  to  be  used.  The  mod-2  adder  outputs  are 
arranged  in  groups  of  2  bits  to  provide  n  4-ary  output  symbols. 
This  program  computes  distances  and  weights  of  each  error  path 
for  constraint  lengths  up  to  K»15  and  input  data  sequence 
lengths  up  to  31  bits. 

The  convolutional  encoder  is  specified  by  selecting  the 
shift  register  length,  the  input  symbol  radix  and  the  output 
symbol  description.  Each  output  symbol  is  determined  from  the 
output  symbol  radix  and  the  bits  in  the  shift  register  with  tap 
connections  which  are  used  in  each  mod-2  summation.  The  tap 
specification  is  used  to  generate  a  mask  which  is  logically 
"and-ed"  with  the  data  and  mod-2  added  to  produce  the 
appropriate  bit  in  the  output.  The  output  distances  are 
accumulated  as  each  data  sequence  is  shifted  through  to  produce 
the  codeword  distance  for  that  sequence.  The  number  of  "ones" 
in  the  data  sequence  is  added  to  the  weight  accumulator  for  the 
corresponding  codeword  distance  to  determine  the  total  weight. 
The  program  uses  a  modified  exhaustive  search  of  all  possible 
patterns  of  the  specified  sequence  length.  Path  errors  begin 
when  the  path  selected  by  the  decoder  differs  from  the  correct 
(all-zero)  path,  so  that  sequences  which  begin  with  all  "zeros" 
will  have  already  been  counted  when  the  coder  used  the  same  data 
sequence  with  the  leading  "zeros"  removed.  For  this  reason, 
these  data  sequences  (which  are  multiples  of  the  shift  register 
radix)  are  omitted.  Another  pattern  can  exist  in  the  data 
sequence  which  does  not  contribute  to  the  result,  namely  any 
data  pattern  which  causes  the  shift  register  to  enter  the  all 
"zeros"  state.  This  can  occur  for  data  sequence  lengths  which 
exceed  the  constraint  length  of  the  decoder,  and  in  essence 
cause  the  decoder  to  leave  the  correct  path  and  return  to  it 
more  than  once.  This  is  not  a  single  path  error  and  these  data 


are  not  included  in  the  weight  calculation. 

After  calculating  the  weights  for  each  code  word  distance 
in  the  selected  sequence  length,  the  program  lists  the  weights 
and  input  data  bit  patterns  associated  with  all  distances  found. 
The  truncated  generating  function  can  then  be  determined 
directly  from  these  data.  Table  A-l  shows  the  results  of  this 
procedure  for  a  K»7,  R-l/2  code  using  a  data  sequence  length  of 
18  bits.  This  code  produces  one  quaternary  output  symbol  for 
each  binary  input  symbol.  The  taps  used  for  each  digit  of  the 
output  symbol  are  listed  in  binary  form  and  represent  an  encoder 
structure  as  shown  in  Figure  4  (in  body  of  report).  From  Table 
A-l  we  see  that  this  encoder  has  weights  7,  39,  104,  and  352  for 
distances  7,  8,  9,  and  10  respectively,  as  shown  in  Table  I  of 
the  report.* 


*  This  code  weight  structure  replaces  a  previously  published 
code  weight  structure  [4]  which  was  found  to  contain  numerical 
errors. 
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SEQUENCE  LENGTH  -  18  BITS 
NUMBER  OF  WEIGHTS  LISTED  -  4 

SHIFT  REGISTER  LENGTH  -  7  BITS,  RADIX  - 

NUMBER  OF  OUTPUT  SYMBOLS  -  1,  RADIX  -  4 

SYMBOL  •  1: 

TAP  »  1  -  0111111 
TAP  #  2  -  1101101 


DISTANCE 

7 


NUMBER  OF  PATHS 
4 


TOTAL  NUMBER  OF  ONES  (WEIGHT) 
7 


BIT  PATTERNS 


DISTANCE 

8 


NUMBER  OF  PATHS 
10 


TOTAL  NUMBER  OF  ONES  (WEIGHT) 
39 


BIT  PATTERNS 

111 

10111 

11001111 


1011 

11011 

101110101 


1111 

100101 


10101 

1001001 


DISTANCE 

9 


NUMBER  OF  PATHS 
19 


TOTAL  NUMBER  OF  ONES  (WEIGHT) 
104 


BIT  PATTERNS 

11111 

1001111 

1100111 

10111111 

1011101011 


101111 

1011011 

1110101 

100100101 

1100111111 


110111 

1011101 

10010101 

110011111 

1101110101 


111111- 

1011111 

10100011 

1001001001 


DISTANCE 

10 


NUMBER  OF  PATHS 
62 


TOTAL  NUMBER  OF  ONES  (WEIGHT) 
352 


BIT  PATTERNS 

1101 

100011 

101101 

1000011 

1100011 

10011011 

10111011 

11110101 

101110111 

1001000011 

1011001111 

10100011011 

11001111111 

100101110101 

1011101011011 

11001111001111 


10011 

loom 

110001 

1001011 

1100101 

10011111 

11011011 

100100011 

101111111 

1001001111 

1111001111 

10111010101 

11011001111 

101101110101 

1011101110101 

101110101110101 


11001 

101001 

110011 

1010001 

1111111 

10110101 

11011101 

100111111 

110001001 

1010001001 

10010010101 

10111010111 

11011101011 

lrooiiiioioi 

10100011001111 


11101 

101011 

110101 

1010101 

10010111 

10110111 

11101011 

101000011 

110011011 

1010100011 

10010100011 

11001110101 

100100100101 

1001001001001 

.0111010100011 


Table  fl-1  Weight  structure  for  K»7  code. 
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APPENDIX  B:  CHERNOFF  BOUND  FOR  BINARY  FSK  ERROR  PROBABILITY  ON 

r'SAYcmH  fading  caam: 


The  optimal  receiver  for  noncoherent  binary  FSK  on  a 
Rayleigh  fading  channel  consists  of  two  pairs  of  matched  fil¬ 
ters,  one  pair  matched  to  each  frequency  tone.  In  Figure  B-l. 
we  see  that  each  matched  filter  pair  consists  of  one  filter 
matched  to  the  in-phase  component  and  one  filter  matched  to  the 
quadrature  component.  The  output  of  each  matched  filter  pair 
is  a  pair  of  i.i.d.  random  variables  with  probability  density 

O 

G(0,<r  ).  For  the  "correct"  filter  pair  (assuming  signal  1  is 

sent)  the  outputs  x^  and  y^  have  variance  =  °  ^  —  .  and  for 

the  "incorrect"  filter  with  outputs  x„  and  v. ,  the  output 
2  N0  r  %  JX 

variance  is  0-2  *  ^  .  (Here,  E  is  the  average  energy  of  the 

received  signal) . 


2 


Fig.  B-l  —  Optimal  receiver  for  noncoherent  BFSK 
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For  the  receiver  configuration  shown  in  Figure  B-l, 
the  output  statistic  is  w,  the  difference  between  the  squares 
of  the  envelopes  of  each  filter  pair  output.  If  signal  1  is 
sent,  an  error  will  occur  if  w>0.  The  resulting  probability 
of  error  Pg  is  Pr  jw>0|  so  that 

oo 

p  *  /  p(w)dw  (B-l) 

•'o 

where  pCw''  is  the  probability  density  function  for  the  output 

random  variable  w.  To  obtain  a  Chernoff  bound  on  P  we  over- 

e 

bound  the  unit  step  u(w)  ■  1  (for  w  >0)  by  the  exponential 
function  e^wJ-oo<w<oo,  where  X  is  a  free  parameter,  (x>0).  This 
is  shown  in  Figure  B-2. 


Fig.  B-2  —  Exponential  overbound  of  unit  step  indicator  function 
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Thus  we  may  bound  the  error  probability  by  a  function  of  A.  , 


r°° 

?e  <  B(X)  =  J  exw  p(w)  dw 


(B-2) 


The  right-hand  side  of  Eq.  (B-2)  is  the  moment  generating  Unction 
for  the  random  variable  w .  The  fact  that  the  error  bound  is  a 
moment  generating  function  is  the  reason  why  the  error  bound 
for  N  uses  of  a  memoryless  channel  is  simply  the  power  of 
the  bound  for  a  single  use.  Tightening  this  bound  by  selecting 
the  free  parameter  x  yields  the  Chemoff  bound  B: 


B  -  BOO  -  "*■>  /  p(„>  dw  -  j*»  (B-3) 


where  the  wiggly  overbar  indicates  the  expectation  operator. 


For  this  particular  channel,  we  may  perform  this  averaging 

over  the  individual  (independent)  components  of 

9  7  2  2  2 

w  »-  x2  +  y2  -  x^  -  y-^  .  The  random  variable  x-^  is  G^.cr^  ), 

so  chat  2 


2  ^2 

°°  -XX,  2a7 


‘l  ,  1  /*  e  1  e2al 

X~TT  J 


iliro^  -oo 


dx^  = 


L+2a12x 


(B-4) 


Similarly 


fl+2<7^  x 


(B-5) 


xx2  xy2 

e  ■  e 


jl-2a2  A 


(B-6) 
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(B-7) 


j^l+2 o^xj  1^1- 2ax  Aj 

[l+(E+H0)x]  [rr^x] 


Letting  AN 


We  find 


A  and  1  + 


B(A) 


<1+0a')(1-X')  l+</3-l)X'-/JAy 


(B-8) 


Differentiating  the  denominator  with  respect  to  a/  and  equating 
to  zero,  we  find 


2(l+g-) 

o 


(B-9) 


Applying  this  to  B(A),  we  obtain  the  Chemoff  bound  for  the 
Rayleigh  fading  channel. 


The  optimal  receiver  for  noncoherent  binary  FSK  on  a 
partial  band  Gaussian  channel  is  the  same  as  shown  in  Figure 
B-l.  For  partial  band  Gaussian  jamming,  there  is  a  probability  p 
of  hopping  into  a  jammed  portion  of  the  band  where  the  noise 
density  is  NQ/p;  otherwise,  the  transmission  is  noise  free. 

During  jamming,  the  (independent)  output  statistics  (for  the 
receiver  in  Figure  B-l)  have  probability  density  functions  as 
follows : 


y2  •  G  (°.  S) 


Following  the  approach  used  in  Appendix  B,  the  Chemoff  bound 
for  the  probability  of  error  (with  worst  case  jamming)  is 

n  _  max  min 
0<p<l  \>0 


/ 

Letting  Nq  -  ~  ,  we  get 


(C-2) 

(C-3) 
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Xx02 


xy22 


Letting  X/  *  XN^  ,  we  find 


e 


(  2  2  2  2\  ~ A 

x2  +  y2  -  *i  -  yi  )  i  i+v 

'  - 77  e 

1-X7  ^ 


-x'  E 


7R; 


Thus 


p  . «  _  max  min 

eS  0<«1  x/>0 


l-x' 


:^L  p£-' 

p  l+x'  ^ 
77  e 


Taking  ^  and  equating  to  zero,  we  find 


i+\'  ^ 


so  that 


P  <  B 
e 


min 

x7  >0 


e'1 _ 1 _ 

^  /a-*') 


Taking  — y  ,  and  equating  to  zero  we  find 

dx 


(C-4) 


(C-5) 


(C-6) 


(C-7) 


(C-8) 


(C-9) 
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Thus,  the  Chernoff  bound  for  the  partial  band  Gaussian  channel  Is 


P  <  B 
e  - 


(C-10) 


This  bound  Is  effective  In  the  range  E/NQ  >  3,  which  Is  the  range  of 
Interest  In  this  application.  For  E/NQ  <  3,  broadband  jamming  (p*l)  is 
the  worst  case  jamming. 


29 


